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4bstract 


We investigated the X-ray scattering signal of highly aligned multilayers of the zwitterionic lipid l,2-dilauroyl-sn-glycero-3- 
shosphatidylcholine containing pores formed by the antimicrobial peptide gramicidin as a function of the peptide/lipid ratio. We 
ire able to obtain information on the structure factor of the pore fluid, which then yields the interaction potential between pores 
n the plane of the bilayers. Aside from a hard core with a radius close to the geometric radius of the pore, we find a repulsive 
ixponential lipid-mediated interaction with a decay length of 2.5 A and an amplitude that decreases with the pore concentration, 
n agreement with the hydrophobic matching hypothesis. In dilute systems, the contact value of this interaction is about SOksT. 
similar results are obtained for gramicidin pores inserted within bilayers formed by the nonionic surfactant pentaethylene glycol 
nonododecyl ether. 


Key words: Gramicidin, Small-angle x-ray scattering. In-plane interaction. Hydrophobic mismatch 
PACS: 61.05.cf, 87.16.dt, 87.15.kt 


1. Introduction 


In the last decades, much effort has been dedicated to the 
understanding of biological membranes, in particular to the 
interaction between membrane proteins and the host lipid 
□ilayer. While the fluid mosaic model [1] described the pro¬ 
teins as free to diffuse in this environment, it was perceived 
very early on [2] that, far from being a neutral background, 
the lipid bilayer can influence protein organization in the 
plane of the membrane and hence many aspects of their ac¬ 
tivity (such as cell signaling and membrane trafficking). In 
particular, cholesterol content was recognized as an impor¬ 
tant variable. For instance, it was put forward as a regulat¬ 
ing parameter for protein partitioning between the plasma 
membrane and the Golgi complex [3]. More generally, the 
homogeneity of cholesterol-containing membranes and the 
biological relevance of the so-called lipid rafts [4] has been 
the subject of extensive research [5]. 

Interestingly enough, the mechanisms involved in 
protein-lipid interaction and in the resulting protein- 
protein interaction mediated by the lipid membrane are 
non-specific, in that no chemical bond is formed between 
a protein and a small number of definite lipid molecules. 
Thus, a detailed understanding and a quantitative char¬ 
acterization of the interaction between lipids and proteins 
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in membranes require that one considers the membrane as 
a many-particle system whose properties are collectively 
determined by the assembly and not only by the chemical 
properties of the individual lipids and proteins [6]. This 
justifies the hope that -notwithstanding the complexity of 
the system- the concepts developed in soft matter physics 
for the understanding of self-assembled systems are op¬ 
erative in this context and that ’simplified’ models can 
yield valuable information. For this reason, a considerable 
body of work dealing with the theoretical modelling and 
numerical simulation of protein-lipid systems appeared in 
the last decades, see [7] for a recent review. These efforts 
are either continuum-elasticity theories or more detailed 
models taking into account the molecular structure of the 
lipid bilayer. 

However, very few experiments attempted to deter¬ 
mine directly the interaction forces between membrane 
inclusions. First among them, freeze-fracture electron mi¬ 
croscopy (FFEM) studies [8-11] yielded the radial distri¬ 
bution function of the inclusions. Comparing the data to 
liquid state theories [12-14] resulted in a hard-core model 
with, in some cases, an additional repulsive or attractive 
interaction. FFEM was not extensively used, undoubtedly 
due to the inherent experimental difficulties; moreover, it 
is not obvious that the distribution measured in the frozen 
sample is identical to that at thermal equilibrium. 

The interaction of membrane inclusions can also be stud- 





ied using small-angle neutron or X-ray scattering from ori¬ 
ented samples, as demonstrated by Huang and collabora¬ 
tors [15-17]. One can thus measure the structure factor of 
the two-dimensional system formed by the inclusions in the 
membrane. Further analysis gives access to the interaction 
potential of the inclusions. In our experience, this method 
is more reliable when several measurements are performed 
along a dilution line (varying density of inclusions) and the 
results are treated simultaneously. As an application, we 
could measure a repulsive interaction both for alamethicin 
pores in DMPC membranes [18] and for inorganic particles 
in synthetic membranes [19]. 

In this paper, we study pores formed by the antimicro¬ 
bial peptide gramicidin D in bilayers with different compo¬ 
sitions. Simultaneously fitting the two-dimensional struc¬ 
ture factor of the pores in the plane of the membrane yields 
the interaction potential between the pores. 

2. Materials and Methods 

2.1. Sample preparation and environment 

The lipid l,2-dilauroyl-sn-glycero-3-phosphatidylcho- 
line (DLPC) was purchased from Avanti Polar Lipids 
Inc. (Birmingham, AL, USA). The antimicrobial pep¬ 
tide gramicidin D and the zwitterionic surfactant N,N- 
dimethyldodecylamine-N-oxide (DDAO) were bought from 
Sigma Aldrich. The nonionic surfactant pentaethylene gly¬ 
col monododecyl ether (C 12 EO 5 ) was bought from Nikko 
Chemical Ltd. (Japan). Without further purification, the 
products were dissolved in isopropanol. The stock solutions 
were then mixed to give the desired molar peptide/lipid 
ratio P/L. The resulting solutions were then dried in vac¬ 
uum and hydrated in excess water (DLPC phases) or up 
to a water content of 20 wt% (DDAO and C 12 EO 5 ). 

The samples were prepared in flat glass capillaries (Vitro- 
Com Inc., Mt. Lks, NJ, USA), 100 /rm thick and 2 mm wide 
by gently sucking in the lamellar phase using a syringe. The 
capillaries were flame-sealed. Areas aligned in homeotropic 
alignment (lamellae parallel to the flat faces of the capil¬ 
lary) formed slowly (over a few months, at room temper¬ 
ature) in DLPC phases. Samples of DDAO and C 12 EO 5 
were aligned by thermal cycling between the lamellar and 
isotropic phases, at cooling rates of about 1 °C/ min. 

2.2. Measurement 

The SAXS measurements were performed at the bend¬ 
ing magnet beamline BM02 (D2AM) of the European Syn¬ 
chrotron Radiation Facility (ESRF, Grenoble, France). The 
photon energy was set at 11 keV. See reference [20] for more 
details. 

The data was acquired using a Peltier-cooled CCD cam¬ 
era (SCX90-1300, from Princeton Instruments Inc., NJ, 
USA) with a resolution of 1340 x 1300 pixels. Data pre¬ 
processing (dark current subtraction, flat field correction. 


radial regrouping and normalization) was performed using 
the bm 2 img software developed at the beamline. 

The incident beam was perpendicular to the flat face of 
the capillary (parallel to the smectic director, which we 
take along the z axis.) Thus, the scattering vector q is 
mostly contained in the (x, y) plane of the layers, and the 
measured scattered signal /(q) probes inhomogeneities of 
the electron density in this plane. Since the bilayers form 
a two-dimensional liquid, the scattering pattern exhibits 
azimuthal symmetry: I = I{q = ]q|). The capillaries were 
scanned in the beam to find well-aligned domains (where 
the intensity of the residual Bragg reflections was as low as 
possible.) 


2.3. Analysis 


The gramicidin pores are dispersed in the lamellar phase 
matrix. Since the “pure” lamellar phase gives a signal con¬ 
fined to the vicinity of the Bragg peaks, from the Babinet 
principle it ensues that the off-axis scattering is the same 
as for a system where the density profile of the lamellar 
phase is subtracted, and one is left with fictitious “pore mi¬ 
nus bilayer” objects in a completely transparent medium. 
Furthermore, as the pores represent a collection of identical 
and similarly oriented objects (up to an azimuthal averag¬ 
ing), the classical separation of the scattering intensity in a 
structure factor multiplied by a form factor can be applied 
[21], yielding: /(q) = S'(q) • E/(q), with: 


S{qz,qr) 


1 

N 


N-1 

exp (-fqrfe) 


( 1 ) 


where N is the number of objects and object “0” is taken as 
the origin of the coordinates. If there is no in-plane order¬ 
ing, S only depends on the absolute value of the in-plane 
scattering vector qr = Jq^ + 'Zy ■ 


The form factor Ff{qr) is the squared modulus of the 
Fourier transform of the electron density p{r) of the scat¬ 
tering object. We computed it from the atomic coordinates 
of the starting structure used in the molecular dynamics 
(MD) simulation of [22], projected onto the z plane. The 
main feature of the density profile is the increased electron 
density at the position of the backbone, which we fitted by 
a radial Gaussian profile, with a peak radius ro = 3.47 A 
and a width w = 1.5 A. The resulting from factor (Fig¬ 
ure 1) is very close to that obtained by [23] from the MD 
simulations of [24]. 

The intensity is divided by the form factor to yield the 
two-dimensional structure factor S{q) of the fluid formed 
by the pores in the plane of the membrane (Figure 2.) The 
uncertainty in the intensity is estimated from the spread of 
the values recorded by the different pixels with the same 
q value. We use standard error propagation through the 
background subtraction and form factor division steps to 
estimate the uncertainty in the structure factor values. 
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3. Results and Discussion 


CT* 


Form factor of the 
gramicidin pore 

O present work 
-Harroun et al. [23] 



Fig. 1. Normalized in-plane form factor of the gramicidin pore (he¬ 
lical dimer configuration) obtained from the atomic configuration of 
deGroot et al. [22] (line and symbols) compared to the form factor 
used by Harroun et al. [23] (solid line). 



Fig. 2. Scattered intensity I{q) (symbols) for a sample with 
PjL = IjT.b and normalized form factor (solid line); see Figure 1. 
Inset: the structure factor obtained by dividing the measured inten¬ 
sity through the form factor. The light symbols in the intensity curve 
(and gaps in the structure factor) correspond to unusable data: at 
small angles (in the shade of the beamstop), at the position of the 
lamellar peaks and at high angles, where a systematic discrepancy 
is observed between the intensity and the form factor. 


The sample alignment is not perfect and sometimes resid¬ 
ual lamellar peaks persist, as seen in Figure 2. The points 
around these positions are discarded (light symbols in the 
I{q) curve and gaps in the S((i) curve.) We also discard 
the points at low g-values (in the shade of the beamstop) 
and those above 0.5 A , where a systematic discrepancy 
between the experimental data and the model form factor 
leads to an artificial oscillation of the structure factor. In 
spite of these limitations, the first peak of the structure 
factor is properly measured for all samples. 


The structure factor curves are plotted as open dots 
in Figure 3 a) for all PjL values, indicated alongside the 
curves. The uncertainty bars are shown for all data points 
(they are generally smaller than the symbol size.) In the fol¬ 
lowing, we will make the simplifying assumption that there 
is no interaction between the pores along z (from one bi¬ 
layer to the next), a result obtained by [17] for gramicidin 
pores in fully hydrated bilayers. This feature is extremely 
important for two reasons: first of all, it allows us to use a 
two-dimensional model, only dealing with the interaction 
within a bilayer; second, it ensures the biological relevance 
of this study. 

3.1. Hard disk model 

The simplest model for the interaction of gramicidin 
pores in membranes is that of hard disks confined in the 
plane. Such an analysis was already performed by [23] for 
gramicidin in DLPC bilayers at P/L = 1/10. As a first 
step, we analysed all the curves using the two-dimensional 
structure factor S'hd(Qr)) given by the simple analytical ex¬ 
pression obtained by [25] (see Eq. 6.8) using the “funda¬ 
mental measure” approach. The details are shown in Fig¬ 
ure 3 for the gramicidin/DLPC system: panel a) displays 
the data and fits, while panels b) and c) show the evolution 
with PjL of the fit parameters, namely the number density 
of the pores ripore and the effective hard-disk radius Rhd- 
Both parameters vary freely during the fit. 

A first observation is that the pore density (symbols in 
Figure 3 b) is in very good agreement with the theoretical 
value calculated using the molar ratio P/L and published 
data for the area per lipid/surfactant molecule and gram¬ 
icidin pore. A similar agreement is obtained for the two 
other systems: gramicidin/Ci 2 E 05 and gramicidin/DDAO 
(data not shown). 

The second fit parameter, the effective hard disk radius 
Rhd is shown as a function of PjL in Eigure 3 c) and as 
a function of the pore density Upore and the area fraction 

b = Hpore X Apore (with Apore = 250 A the area of a gram¬ 
icidin pore) in Eigure 4. Eor comparison, we also plot the 
data point of [23] (open square). From the analysis of the 
structure factor for Gram/DLPC at P/L = 1/10, these au¬ 
thors obtain a hard-disk radius of 13.4 A, which they in¬ 
terpret as the geometric radius of the pore plus one lipid 
layer. This value is somewhat higher, but still compatible 
with our data ^. 

In Figure 4, we also show the fit results for the two other 
systems. For all three membrane compositions, the effec¬ 
tive radius decreases as a function of Upore (although this 

^ This discrepancy might stem from the slightly different way of 
obtaining the interaction radius. Harroun et al. [23] determine the 
most probable nearest-neighbor separation between two pores, al¬ 
ways larger than twice the hard-disk radius, which is the minimum 
separation. 
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Fig. 3. (a) Experimental structure factors (symbols) and hard-disk fits (solid lines) for the gramicidin/DLPC system, at different peptide/lipid 
(P/L) molar concentrations indicated alongside the curves. The theoretical surface fraction of gramicidin pores, 77, is also given, (b) Pore 
density obtained from the fits in (a) (symbols) and theoretically expected value (dotted line), (c) Effective values for the hard disk radius 
obtained from the fits in (a) (symbols) and geometrical pore radius (dotted line). The variation in Rhd with the pore concentration is a 
sign of repulsive interaction, as discussed below. 


decrease is less marked for gramicidin/DDAO.) Clearly, the 
hard disk model is not satisfactory: although it fits very 
well the individual curves, the interaction radius decreases 
as P/L (and hence Upore) increase; as already discussed for 
the case of alamethicin pores [18], this is a sign of an ad¬ 
ditional “soft” repulsive interaction. Briefly, this effect can 
be understood as follows: at low density, this interaction is 
enough to keep the particles far from each other. As the 
density increases so does the pressure, and the particles are 
eventually pushed closer and closer to each other. At high 
concentration, the effective radius should saturate at the 
(impenetrable) “true” core value. This seems indeed to be 
the case for the data in Figure 4, where at increasing den¬ 
sity the curves approach the geometric core radius of 9 A. 

3.2. Additional interaction 


To quantify this tendency, we calculated the structure 
factors for a hard core with radius Rhd = 9A and an 
additional “soft” potential: 


V{r) = u exp 


1 

2 



r > 2 Rhd 


( 2 ) 


where r denotes the distance between the pore centres. 

The structure factor S{q) is now a function of four pa¬ 
rameters: the hard core radius Rhd, the number density 
Upore, as well as the amplitude u and the decay length ^ 
of the additional component. We calculate S{q) using the 
method of Lado [26, 27], implemented as an iGOR Pro 
function ^ . Briefly, the method provides an iterative solu¬ 
tion to the Ornstein-Zernicke equation with the Percus- 
Yevick closure. The accuracy of the procedure was tested 
by comparing the results for purely hard-core systems to 
the analytical formula of Rosenfeld [25]. 

The effect of the additional potential on S{q) is illus¬ 
trated in Figure 5, displaying the position (jmax of the first 
maximum in the structure factor as a function of pore den¬ 
sity for the Gram/DLPC system, as well as for simulated 
structure factors with a hard core Rhd = 9 A (the geomet¬ 
ric core radius of the pore) and an additional component 
(Equation 2) with a decay length ^ = 15 A and various am¬ 
plitudes u (positive and negative). 

The first observation is that, for a purely hard-core inter¬ 
action, the position of the maximum in S{q) varies very lit- 


^ The source code is available from the author upon request. 
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Fig. 4. Effective hard disk radius of the gramicidin pore as a function 
of the density in bilayers with three different compositions. The lines 
are just guides for the eye. 
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Fig. 5. Position of the first maximum ^max in the structure factor 
S{q). Experimental data (A and thick line) for the Gram/DLPC 
system and simulated data (various symbols and lines) for different 
amplitudes of the additional interaction u. For all simulations, the 
hard core radius Rhd = OA and the decay length ^ = 15A. 

tie with the concentration. Adding an attractive potential 
u < 0 only shifts the peak to higher q values (the pores are 


on the average closer to their neighbours) without chang¬ 
ing its slope. When the potential is repulsive m > 0, on the 
other hand, the maximum shifts to lower q values and its 
density dependence becomes steeper, although not enough 
to describe the variation of the experimental data, even 
when the parameters u and ^ are allowed to vary within 
reasonable limits. A preliminary conclusion would thus be 
that the interaction potential contains a “soft” repulsive 
contribution, induced by the membrane, which varies with 
the density. 

To check this conclusion, we made extensive attempts 
to fit the data to a hard-core model plus an additional in¬ 
teraction (Equation 2), where Rhd = 9A is the same for 
all curves, and ripore is fixed at the theoretically expected 
value, given by the dotted line in Figure 3, while u and 
^ are allowed to vary freely. Alternative functional forms 
for the additional interaction (exponential and linear) were 
also tested. None of these attempts yielded satisfactory re¬ 
sults; in particular, parameter combinations that capture 
the density dependence of the peak position q^ax also give 
a marked variation in its width and amplitude, in contrast 
with the experimental data in Figure 3 a), where the shape 
of the peak changes very little as P/L varies by a factor of 
five. 

3.3. Possible origins for the interaction 

To see why and in what manner the membrane-mediated 
interaction should vary with the pore density, we must first 
consider the effect of inclusions on the membrane. Two such 
effects are relevant in this context. The first one is related to 
changes in the membrane thickness and occurs on a “meso¬ 
scopic” scale, affording a continuum elasticity treatment. 
The second effect concerns the way in which a membrane 
inclusion perturbs the configuration of the lipid chains and 
requires a more involved, microscopic description. 

3.3.1. Hydrophobic matching 

One of the simplest (and surprisingly successful) con¬ 
cepts used to interpret protein-membrane interaction is 
that of hydrophobic matching [28, 29], stating that pro¬ 
teins with a certain hydrophobic length (defined by their 
transmembrane domain) are targeted to membranes with 
a matching hydrophobic thickness. When there is a differ¬ 
ence in length between the hydrophobic part h of the pro¬ 
tein or peptide and that of the host membrane (referred to 
as the “hydrophobic mismatch”), the bilayer deforms and 
adapts to the protein (which is generally much more rigid). 
This deformation has a certain lateral extension (in the 
plane of the membrane) and therefore induces an interac¬ 
tion between inclusions when the latter are closer than this 
distance. 

The hydrophobic length of gramicidin was estimated at 
ho = “2,2 A [30]. For DLPC, h = 20.8 A in the pure mem¬ 
brane and h = 22.1 A he for a gramicidin content 
PjL = 0.1 [23], confirming that thinner membranes are 















stretched by gramicidin. The same effect is observed in 
the gramicidin/DDAO system, where NMR measurements 
yield a hydrophobic thickness of 18.4 A for the pure bilayer, 
increasing to 19.4 A for the more concentrated samples 
{P/L = 0.25) [31]. For C 12 E 5 , we estimate the hydrophobic 
thickness in the absence of gramicidin at hcijEs = 18.8 A; 
we have no data on the variation of hciaEs with the pore 
density. 

To first order (and ignoring the effect of the spontaneous 
curvature of the monolayer, which can have non-trivial ef¬ 
fects [32, 33]), the elastic energy is expected to scale as 
{h — hcY [34]. It will therefore decrease as the gramicidin 
concentration increases, and vanish when the hydrophobic 
thickness of the membrane equals that of the bilayer. 

3.3.2. Changes in lipid ordering 

On the other hand, even in the absence of hydrophobic 
mismatch, the presence of an inclusion restricts the config¬ 
uration of lipid chains in its vicinity [35-38]. For instance, 
Lagiie et al. [37, 39] used the lateral density-density re¬ 
sponse function of the hydrocarbon chains obtained from 
MD simulations of pure bilayers to determine the interac¬ 
tion between “smooth” hard cylinders embedded in the bi¬ 
layer. They considered three values of the cylinder radius, 
up to 9 A, which is precisely the geometric radius of the 
gramicidin pore. For the largest radius, the long-range in¬ 
teraction is repulsive for all the investigated lipids (DMPC, 
DPPC, POPC, DOPC), with an additional short-range at¬ 
traction in the case of DMPC. However, no specific predic¬ 
tions are available for DLPC. Furthermore, it is not clear 
how the interaction varies as a function of inclusion con¬ 
centration. 

3.4. Complete model 

In the following, we adopt the hydrophobic matching 
model and assume that the mismatch {h — he) varies lin¬ 
early with the in-plane concentration of pores Upore over 
the investigated range (1/25 < PjL < 1/5, or rimin = 
0.55 < Upore < ?T-max = 1-7 X 10“^ A ) and vanishes at 
the highest value. We model this elastic interaction by an 
exponential ^ 

T/i' ^ TT (~ ?^pore\ ( T ~ 2R\ 

V{r)=Uo[ -^ exp--— (3) 

\ ^max ^min / \ S / 

where the free parameters are the prefactor Uq, which is 
also the interaction amplitude at the lowest concentration 
nmin, and the decay length taken as independent of the 
concentration. No effort was made to include three-body 
interactions. Since the additional potential vanishes at the 
highest pore concentration, in the following we take as hard 
core radius the best fit at this value, namely 9.5 A, rather 


® We also performed tests using different functional forms for the 
interaction potential, but the fits are less satisfactory (see the ap¬ 
pendix). 
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Fig. 7. Experimental structure factors (symbols) and fits (solid 
lines) with a hard-disk plus the repulsive interaction (3) for the 
gramicidin/Ci2E05 system. The pore density is indicated alongside 
each curve (the four curves correspond to the solid squares in Figure 

4 ). 

than the geometric value of 9.0 A. This adjustment is not 
very significant, but is required for coherence of the model 
and it ensures that the additional component is not over¬ 
estimated. 

For the Gram/DLPC system, the best fit with this model 
is obtained for Uq = 31.5 ± 10 ksT and ^ = 2.5 ± 0.5 A. 
The details are given in Figure 6. See the Appendix A for 
details on the error estimate. 

A similar result is obtained for the Gram/Ci 2 E 05 sys¬ 
tem: C/q = 27± 10 ksT and ^ = 2.75 ±0.5 A (see Figure 7). 
As above, the mismatch {h — he) is taken to vary linearly 
from a maximum at the lowest pore concentration to zero 
in the most concentrated sample. 

For the Gram/DDAO system, fitting attempts point to¬ 
wards a much lower membrane-mediated interaction (both 
in amplitude and decay length), coherent with the low vari¬ 
ation in effective hard disk radius (open dots in Figure 4.) 
Clearly, the interaction of the pores is closer to a pure hard¬ 
disk model than for the other bilayer compositions. How¬ 
ever, a quantitative analysis is difficult in view of the low 
quality of the fits. Further experimental data is needed on 
this system. 

4. Conclusion 

We showed that the interaction between gramicidin pores 
inserted within DLPC bilayers can be well described by a 
hard disk potential, with a range close to the geometric di- 










Fig. 6. (a) Experimental structure factors (symbols) and fits (solid lines) with a hard-disk plus the repulsive interaction (3) for the grami- 
cidin/DLPC system (data as in Figure 3). (b) Top: “Effective fraction” of the interaction amplitude as a function of npore: fr = ( \ ^ 

see Equation (3). Bottom: Position of the structure factor maximum gmax vs. ripore, for the experimental data (symbols) and for the 
fits (solid line), (c) Interaction potential V(r) used for the fits in panel (b). The amplitude corresponds to the lowest value of ripore, 
«min = 0.55 X 10-3 


ameter of the molecule and an additional repulsive interac¬ 
tion, whose amplitude decreases with increasing pore con¬ 
centration. A similar result is obtained for pores inserted 
within C 12 EO 5 bilayers. 

The decrease of the interaction amplitude with pore con¬ 
centration is consistent with the hydrophobic matching hy¬ 
pothesis, whereby the increase in bilayer thickness (well- 
documented for Gram/DLPC systems [23]) helps “accom¬ 
modate” the peptide and leads to a decreased interaction. 
The decay length of the interaction is found to be of the or¬ 
der of 2.5 A, well below the value predicted by continuum 
elastic models [32, 40] and close to the characteristic length 
given by more microscopic models [41]. However, consider¬ 
ing that the amplitude of the interaction is quite high, the 
interaction range -defined intuitively as the distance over 
which the pores ‘see’ each other- does extend over several 
decay lengths; see Appendix B for a more detailed compar¬ 
ison. 

An intriguing consequence of the interaction decreasing 
with the concentration of pores is that the structure factors 
of the two-dimensional fluid that they form in the plane of 
the membrane changes relatively little with the area density 


(as compared to e.g. purely hard-core interactions). One 
could say that the degree of “liquid order” remains sensibly 
the same. It is tempting to speculate on the biological sig¬ 
nificance of this feature. For instance, does the activity of 
gramicidin require a certain degree of correlation between 
pores? We hope that systematic studies of the interaction 
in various conditions will help answer this question. 


Appendix A. Data fitting 

The first observation is that the (statistical) uncertainty 
in the structure factor S{q)^ determined by radial regroup¬ 
ing, background subtracting and then dividing by the cal¬ 
culated form factor, is much too small. Clearly, the dis¬ 
crepancy between the fit functions and the data is mostly 
due to systematic effects, which can be related to errors 
in determining S{q) from the measured intensity or to the 
inadequacy of the fit functions. We therefore used a more 
realistic uncertainty value, cr = 0.1 taken as a constant for 
all data points. The goodness-of-fit function is then of 
the order of 2 (per data point.) 
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Fig. A.l. The best results for the interaction potential V{r) within 
each model class. The respective values are also indicated. 

Several trial functions for the interaction potential V (r) 
were tested: exponential decrease, gaussian function and 
linear slope. The exponential model presented in Figure 6 
yields the best fit. The best potentials obtained for each 
class of model are displayed in Figure A.l. The respective 
values are indicated alongside the curves. 

Once the fit is obtained, the confidence intervals should 
be determined for each parameter. For the exponential 
model, which yields the best fit, we computed X^{Uo,0 
for a wide range of parameters, to make sure that we have 
indeed found the global minimum and to visualise the de¬ 
pendence of on the parameters. This graph is plotted in 
Figure A.2. Clearly, ^ is within the range 2.5±0.5 A. On the 
other hand, the interaction amplitude is less well defined; 
we take the conservative estimate Uq = 31.5 ± 10 ksT. 

Appendix B. Comparison of the interaction range 

As mentioned in the Conclusion, ^ is not entirely repre¬ 
sentative, since the interaction amplitude can remain signif¬ 
icant over several decay lengths. This is an obstacle to com¬ 
paring the different theoretical models between them and 
with the experimental data, unless they can be described 
by a common functional form. Another -less precise- op¬ 
tion is to dehne an arbitrary cutoff amplitude and to de- 
hne the interaction range as the separation beyond which 
the interaction falls below this cutoff. We did this for our 
model with two cutoff values, Uc = 0.5 and 1 /cbT, respec¬ 
tively. The resulting cutoff ranges are shown in Figure 
B.l (symbols and dotted lines). For reference, the decay 
length ^ is also plotted as solid line. 

In continuum elasticity models, the characteristic length 
is A = [40], with K and B the bending and com¬ 

pression moduli of the bilayer, respectively. These authors 
conclude that a characteristic length in the range 8.5 < 
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Fig. A.2. Goodness-of-fit function as a function of the parameters 
Uq and 


A < 12.5 A (shown in Figure B.l as gray shading) accounts 
for the behaviour of gramicidin in DMPC bilayers. We also 
plotted as dashed line the characteristic length used by 
Aranda-Espinoza et al. [32] (corresponding to /3 = 10 in 
their choice of parameters). Clearly, these values are much 
larger than our value for the decay length, ^ = 2.5 A, with 
A/^ between 3 and 5. 

In order to compare the effective interaction range, one 
would need to estimate the cutoff ranges for the theoretical 
results. Although this data is not available in the publica¬ 
tions, it appears that the interaction extends at least over 
the thickness of one monolayer (shown in Figure B.l as solid 
line), which is also quoted as the dominant scale length in 
reference [42]. This length scale is clearly larger than our 
estimates for the cutoff range, although the discrepancy is 
less marked than between A and 
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Fig. B.l. Distance scales relevant for the interaction: the decay length 
^ measured in this paper and the monolayer thickness (both shown as 
solid lines), the cutoff distances Vc for two cutoff amplitudes, 0.5 and 
1 k^T (symbols and dotted lines) and the values of the characteristic 
length A as given by [40] (gray shaded range) and estimated from 
[32] (dashed line). 
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